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Jiirg Frohlich and Tom Spencer have profoundly shaped the interface between physics 
CS| ■ and mathematics. It is therefore a particular pleasure to dedicate to them this article on 

^ ! a chapter of mathematical physics which is far from being closed. 
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t~^ ■ Abstract. The derivation of the Nordheim-Boltzmann transport equation for weakly in- 

^ . teracting quantum fluids is a longstanding problem in mathematical physics. Inspired by 

(^ ! the method developed to handle classical dilute gases, a conventional approach is the use 

of the BBGKY hierarchy for the time-dependent reduced density matrices. In contrast, 
our contribution is motivated by the kinetic theory of the weakly nonlinear Schrodinger 
equation. The main observation is that the results obtained in the latter context carry over 
directly to weakly interacting quantum fluids provided one does not insist on normal or- 
der in the Duhamel expansion. We discuss the term by term convergence of the expansion 
and the equilibrium time correlation (a(t)*a(0)). 
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1 Introduction 

With the discovery of quantum mechanics, evidently, Boltzmann's kinetic theory of rari- 
fied gases had to be revised. The modification turned out to be minimal, in a certain sense. 
Only the classical differential scattering cross section had to be replaced by its quantum 
version, which thus depends on whether the gas particles are fermions or bosons. Other- 
wise the structure of the equation remains unaltered. In particular, the stationary solutions 
are still the Maxwellians. In addition, a natural task was to investigate the kinetic regime 
for weakly interacting quantum fluids with no particular restriction on the density. Time- 
dependent perturbation theory, or even the more basic Fermi's golden rule, provides a 
convenient tool. Since the interaction is weak, the differential cross section appears now 
only in the Born approximation. The resulting quantum kinetic equation has a cubic colli- 
sion operator, rather than a quadratic one as in Boltzmann's work. As a consequence, the 
Maxwellians exp[— /5(cl'(/c) — /i)] are no longer stationary. The correct stationary solutions 
are of the form (exp[P(uj{k) — /x)] =1= 1)^^ in accordance with the statistics of the quantum 
fluid under consideration. 

The quantum transport equation was first written down by Nordheim 01 in a paper 
submitted on May 30, 1928 and published a few weeks later. His work is an ingenious 
guess, supported by an H-theorem and by the physically expected stationary momentum 
distributions. Later on more systematic derivations followed [12|. In the kinetic literature 
the transport equation mostly carries the names of Uehling and Uhlenbeck [3J. In their 
1933 paper they study properties of the transport equation, in particular its linearization 
around equilibrium and the long time hydrodynamic approximation. 

There is a fundamental difference between weakly interacting classical and quantum 
fluids in the kinetic regime. For two classical particles a small interaction results in a 
small change of the relative momentum. Therefore the collision operator is not an inte- 
gral operator, as in Boltzmann's classic work, but it is a nonlinear differential operator, 
apparently first realized in 1936 by Landau [4J in the context of plasmas and fluids with 
long range interactions. On the other hand, two weakly interacting quantum particles will 
pass through each other with probability 1 — (9(A^), A the coupling constant, and will 
s-wave scatter with probability (9(A^). Therefore the collision operator is still an integral 
operator. 

In our notes we will consider only quantum fluids interacting through a weak pair 
potential. 

The conventional formal derivation of the transport equation proceeds via Fermi's 
golden rule. Over the years there have been many attempts to improve on the argument. 
On the theoretical side we mention in particular van Hove [5], Prigogine (6\, and Hugen- 
holtz [7|. With the work of Lanford [8] on the microscopic justification of the classical 
Boltzmann equation and the work of Davies |l9l on the weak coupling limit for a small 
quantum system coupled to a large heat bath, the derivation of kinetic equations was rec- 
ognized as a problem of interest to mathematical physics. In fact, it is already stated 
verbally in Hilbert's famous collection of problems as problem No. 6 [TTOll . The derivation 
of quantum kinetic equations still stands as a challenge. We will explain its current status 



in due course. 

Regarding a quantum fluid as starting point, there are two in essence orthogonal ways 
to proceed with a semiclassical approximation. Within the particle picture, it is natural to 
take the limit of a large mass which leads to classical point particles interacting through 
a pair potential. On the other side, from the point of view of operator-valued fields, the 
natural limit is a weak, long range (on the scale of a typical inter-particle distance) poten- 
tial which leads to the nonlinear Schrodinger equation, also known as Hartree equation. 
Both limits can be formulated as an Egorov theorem for operators [[m[T2l[T3l[T4l . In the 
former case one is back to the model studied by Boltzmann while in the latter case the 
kinetic issue concerns a weakly nonlinear wave equation as already studied by Peierls |l2l. 

In filSil we investigate the weakly nonlinear Schrodinger equation and develop a ma- 
chinery for dealing with the high-dimensional oscillatory integrals as they arise in the 
Duhamel expansion of the solution of that equation. The goal of this paper is to explain 
how these novel techniques can be used for the derivation of quantum kinetic equations. 
In fact, we will provide mostly sketches of the argument and prove only a few crucial 
points. Otherwise we would be overwhelmed by technical issues. But we do clearly 
indicate where, to our understanding, there are gaps and new ideas will be required. 

To give a brief outline: Weakly interacting lattice bosons and lattice fermions are in- 
troduced in Section[2]together with the Duhamel expansion of the time evolution operator. 
In Section [3] we explain how this expansion differs from the expansion resulting from the 
quantum BBGKY hierarchy and in Section |4] we make a comparison with wave turbu- 
lence for the weakly nonlinear Schrodinger equation. The main body of novel results are 
in Sections |5] to |7] where we discuss equilibrium time correlations and spatially homo- 
geneous nonequilibrium states in the kinetic limit. In an appendix we summarize a few 
properties of spatially homogeneous quantum kinetic equations. 

Acknowledgments. We would like to thank Laszlo Erdos for many illuminating discus- 
sions on the subject. The research was supported by a DFG grant and by the Academy of 
Finland. The work was partially done at the the Erwin Schrodinger Institute for Mathe- 
matical Physics (ESI), Vienna, Austria, whom we thank for their hospitality. 



2 Weakly interacting quantum fluids, Duhamel expan- 
sion 

We consider quantum particles on the (i-dimensional lattice %^ with lattice points x E Z'^. 
Later on there will be conditions which require d > 3. The particles "hop" independently 
to nearby sites. The hopping amplitude for the relative displacement is given by a : 
Z'^ ^ M. with a{x) = a{—x). a is of compact support, i.e., a{x) = for \x\ > R with 
suitable R. The Fourier transform of a is the dispersion relation uj(k) = a{k). Clearly 
uj(k) = uj{k)* = uj(—k) and to is real analytic. This by itself will not be enough and 
further conditions, originating from the analysis in [jTSll . will be imposed. The particles 
interact through the weak pair potential W, V : Z,'^ ^ 'K, V{x) = V{—x). V is assumed 



to be of compact support and A > 0, A ^ 1. Particles are either bosons or fermions. We 
introduce the corresponding annihilation/creation operators a(x), a{x)*, x G Z'^. In the 
case of bosons they satisfy the commutation relations 

[a{x),a{y)*]=6^y, [a{x) , a{y)] = , [a(a;)*, a(y)*] = , (2.1) 

while in the case of fermions they satisfy the anticommutation relations 

{a{x),a{yy} = 6^y, {a{x) , a{y)} = , {a(x)*, a(?/)*} = , (2.2) 

x,y e Z'^. Here [A, B] = AB - BA and {A, B} = AB + BA. With this notation the 
Hamiltonian of our system of particles reads 

H = y^ a{x — y)a{x)*a{y) + ^A >^ V{x — y)a{x)* a{y)* a{y)a{x) . (2.3) 



if at A = is quadratic and denoted by i/har- Clearly the number of particles, N, is 
conserved, 

N=J2 a{xya{x) and [H, N] = . (2.4) 

For / : Z'^ — > C we introduce its Fourier transform as 

nk)=J2f(^>~'''''''- (2-5) 

with k E T'^ = [— |, ^Y, the rf-dimensional unit torus. The physical momentum is 27rk, 
but our convention has the advantage of minimizing the number of 27r prefactors. The 
inverse transform to (|2.5I) reads 

f{x)= [ dkf{ky'-'-^. (2.6) 

Sometimes it is more convenient to have a function, say W , defined on T'^. Then its 
inverse Fourier transform will be denoted by W . With this notation we have in momentum 
space 

H= dkuj{k)a{k)*a{k) -\ — A / dkidk2dksdk4S (ki + k2 — k^ ~ k^) 

Jjd 2 J(jd-^4 

xV{k2 - k3)d{kiyd{k2yd{ks)d{ki) . (2.7) 

Remark 2.1 (Infinite volume limit). The Hamiltonians (12.31) . resp. (|2.7I) . define a unitary 
dynamics for initial states which have a bounded number of particles. To be more precise, 
we introduce the Fock space 

oo 
^=05n, ;^n = ^2(Z")®r (2.8) 

n=0 



and the projections Pg, 9 = ±. P^^ is the subspace of symmetric and P^^ is the subspace 
of antisymmetric wave functions. Since [H, N] = 0, Pg'^n is invariant and H restricted 
to Pg'^n is a self-adjoint and bounded operator for each n. Of particular interest will 
be translation invariant initial states. Then the number of particles is infinite and the 
expression (12.31 ) is only formal. Thus one first has to restrict the fluid to a box A with 
periodic boundary conditions. In the kinetic limit its side-length must be of order A^. 
Thus, for constant density, the average number of particles is finite and of order A^^"'. 
Hence, for fixed A the dynamics is well defined. Conceptually one would prefer to first 
let A t Z"' and then A ^ together with the appropriate rescaling of space-time. This 
program can be carried through for fermions, in which case the operators a{x), x G 
Z'^, are bounded and the dynamics is defined as a group of automorphisms on the C*- 
algebra of quasi-local observables [fT6l|. Thus one can work directly at infinite volume. 
For lattice bosons the dynamics at infinite volume is not so well understood. For our 
purposes the infinite volume limit is not a central issue. One can first write down the 
Duhamel expansion in finite volume A and then establish that, for the particular initial 
state and the observables of interest, the estimates are uniform in A, hence hold when 

A T z*^. 

Remark 2.2 (Stability). We will study equilibrium time correlations in the kinetic limit. 
For this the /3-KMS state (■)/3,a for H — fiN, fi the chemical potential, has to be well- 
defined and, hence, V has to be a thermodynamically stable potential. For lattice fermions 
no further conditions are needed. For lattice bosons a simple sufficient condition would 
be 

y > . (2.9) 

For a discussion of equilibrium states at small A we refer to IITTlfTSl . In the kinetic limit 
thermodynamic stability does seem to play a role. <0 

Remark 2.3 (Continuum limit). In IfTSl we study a classical field theory by regarding 
a(a;), a{x)* as commutative complex-valued field. To avoid ultraviolet divergence, it was 
necessary to introduce a spatial discretization, as will be explained in more detail in Sec- 
tion HI Since we plan to transcribe the results from flSJ to quantum fluids, we stick to 
the lattice theory. Thus our model is appropriate, for example, for a fluid of electrons in a 
crystal background potential and for bosons in an optical lattice. To understand whether 
a continuum limit is feasible, we recall that 

{aikraik'))(,,o = 6ik-k')W^ik), W^ik) = ^^^^^^^^^ ^ ^ . (2.10) 

Here /5 is the inverse temperature, P > 0, and /i is the chemical potential chosen such that 
uj(^k) — ji > (and suppressed in our notation). We adopt the standard convention WL 
where 9 = 1 stands for bosons and 6* = —1 for fermions. 

In the continuum limit, i.e., for a fluid in empty space, one would replace the position 
space Z'^ by R"' and a* by —A, implying the dispersion relation uj{k) = k"^. Clearly in 
this case 

f Wf{k)dk <oo, (2.11) 



indicating that there is no ultraviolet divergence. Indeed, as discussed in ifTTl [TSl . the 
equilibrium state for small A, and with restrictions on the density, is well-defined. It would 
be of considerable interest to find out whether the analysis in [|T5]| could be extended to 
weakly interacting quantum fluids in M'^. (} 



iHt A^-iHt 



A{0) = A. Using ^A{t) = i[H, A{t)], one obtains the evolution equation 



For an arbitrary operator A we define its Heisenberg evolution by A{t) = e Ae 

dt 

di 



a{x, t) = -i ^ a{x - y)a{y, t) - iA ^ V{x - y)a{y, tya{y, t)a{x, t) , (2.12) 



which in momentum space becomes 

d 



d{ki,t) = — icij(fci)a(/ci, t) — iA / dk2dk-idki5{ki + k2 — k^ — k^) 

xV{k2 - k3)d{k2,tyd{k3,t)d{ki,t) . (2.13) 



Equations (12.121) and (12.131) will be solved as a Cauchy problem given the initial, time 
t = 0, state. In general one should allow for spatial variation on the scale A^^. To keep 
matters simple we will restrict ourselves to a spatially homogeneous situation. Thus the 
initial state (•) is assumed to be translation invariant with good clustering properties. In 
particular 

{d{kyd{k')) = 5{k - k')W{k) (2.14) 

for some smooth W . In addition, (■) has to be gauge invariant, i.e., 

{e'^^Ae-'^^) = {A) (2.15) 

for ?9 G [0, 2ii] . One quantity of interest will be the two-point function at time t. Since the 
dynamics preserves translation and gauge invariance, it must necessarily be of the form 

{d{k, tyd{k', ty = 6{k - k')Wx{k, t) , (2.16) 

Wx > 0, Wxik, 0) = W{k), together with {d{k, t)) = 0, {d{k, t)d{k', ty = 0. 

Kinetic theory studies Wx{k,t/\'^) for small A. A special role will be played by 
quasifree states, which we define first. 

Definition 2.4 (quasifree state). Let pi be a self-adjoint operator on i,2{'^'^) with integral 
kernel pi{x,y). For bosons we impose pi > while for fermions < pi < 1. A 
gauge invariant state {■) is called quasifree with correlator pi, if for all Xi,yj G Z"^, 
i = I, ..., niyj = 1, ..., n, it holds 

m n 

((JJa(xi))*(JJa(yj))) = 5mndete{pi{xi,yjyj..^^^^^. (2.17) 

i=i i=i 



Here the operators are ordered from left to right as they appear in the Yl symbol and 
the empty product is interpreted as 1. det_ is the usual determinant, while det+ is the 
permanent. Translation invariance is reflected by pi(x, y) depending only onx — y, resp. 
in momentum space by (I2.14I ). 

Kinetic theory puts forward a rather simple picture of the dynamics for small A: On the 
microscopic time scale the initial W{k) does not change while the unperturbed dynamics 
forces the state to become quasifree. To say, after a time t which is short on the kinetic 
scale and long on the microscopic scale, for arbitrary n E N, the higher moments are 
approximately of the form, W denoting the inverse Fourier transform, 

n n 

{{lla{x,,t)y{lla{y,,t)))^dete{W{x.-y,))^^^^^_^ (2.18) 

i=i j=i 

with all other moments vanishing. On the kinetic time scale, t = (9(A^^), Wx{k,t) 
changes while preserving the quasifree property. Of course, only in the limit A ^ we 
obtain such a strict separation of the two space-time scales. The initial time slip with 
the dynamics generated by //har was studied by Ho and Landau [fT9ll . Under suitable 
assumptions on cu, they prove quasifreeness in the limit t -^ oo provided the initial state 
is £i-clustering. 

We are mostly interested in the kinetic time scale and thus impose the state (•) to be 
quasifree to begin with. As argued in EOl . see also [[2T]|. if one assumes, up to small 
errors, the state still to be quasifree at the long time t = A^^r, r = 0{1), then it is not 
too difficult to determine the (approximate) evolution equation for Wx{k, t). 

To study {a(k, t)*a(k', t)) the only method currently available is to expand the expec- 
tation value of interest with respect to A, which is achieved through the Duhamel expan- 
sion. It will be convenient to work in the interaction representation and we introduce 

a(fc,l,t) = e''^(^)*a(fc,t), a(fc, -1, t) = e-^'^('=)*a(-fc, t)* . (2.19) 

Then (12.131 ) becomes 



d 
di 



d(ki,a,t) = —iXa / dk2dk'idki5{ki — k2 — k^ 



xi((l + a)Vik2 + ks) + (1 - a)Vik, + k,)) 

X exp [ — it( — auj{ki) — uj{k2) + cruj^ks) + u;(fc4))] 

xd{k2, — 1, t)a(/c3, a,t)d{k4^, 1, t) . (2.20) 



For a product it holds that 



1 n n in~l , n 

-T:Yld{kj,aj,t) = ^{Yla{kj,aj,t))—d{k^,a^,t){ JJ d{kj,aj,t)) , (2.21) 

j=l m=l j=l j=m+l 
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Figure 1: A Feynman diagram with no = 2, n = 3. Up (down) arrows correspond to 
parity +1 (—1). The operator ordering at t = is a*a*aa*aa*aa. The interaction history 
is (1,4,1). 

which, using (|2.20l) . integrates in time to 

n n n „i „ 

1 I a(fcj, o"j,t) = 11 a(fcj, cTj) — iA ^^ (Jm ds dk'2dk'^dk'^ 

i=l i=l m=l J^ h-^"? 

x5{k^ -k'^-k'^- A;^)|((l + am)V{k'2 + k'^) + (1 - a^)V{k'^ + k'^)) 

X exp [ — is{—am^{km) — ^{k'^) + amUj{k'^ + ct;(A;4))] (2.22) 

m— 1 n 

x( JJ^ a(A;j,(Tj,s))a(A;2,-l,s)a(A;3,(Tm,s)a(A;4,l,s)( JJ^ a(A;j, cij, s)) . 

j=l j=m+l 

Here the products are ordered from left to right with increasing label j. 

Iteration of (12.221 ) yields the Duhamel expansion for products of the form 



no 



1 I a(fcj, (Jj,t) , t > . 



(2.23) 



The expansion is most concisely organized through Feynman diagrams, compare with 
Fig.[Il We explain a generic term. 

We denote the final time by t, t > 0. The initial time is 0. There will be n "collisions" 
or interactions, n > 0, which in the Feynman diagram are represented as fusions. We 
subdivide [0, t] into n + 1 time slices with index i E /o,n = {0, . . . , n}. Omitting the zero 
we set /„ = {1, . . . , n}. The i-th slice has length Sj. Thus ^"=o ^i = t and the slices are 
[0, So], . . . , [so + . . . + s„_i, t]. To describe the fusions we start in slice with hq + 2n 
line segments, see Fig. [T] for the orientation of the time axis. At time t = sq exactly one 
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triplet of neighboring line segments fuses into a single line segment, while all remaining 
line segments are continued vertically. Thus in time slice 1 one has exactly uq + 2n — 2 
line segments. At time t = sq + si exactly one triplet of neighboring line segments fuses 
into a single line segment, etc.. The last fusion is at t = sq + . . . + s„_i and in slice 
n one has tiq line segments. Sometimes it is convenient to read the Feynman diagram 
backwards in time. The only change is that "fusion" is turned into "branching". We now 
label the line segments by momenta fcj ., G T'' and by parities aij E { — 1,1}. The first 
index is the label of the slice. In the i-th slice we label the momenta and parities from 
left to right by j = 1, . . . , mj with rrii = uq + 2n — 2i. The corresponding index set is 
denoted by Tn-no = {{hj)\0 '^i<n,\<j< nii}. An interaction history is denoted by 
(, = {£i, . . . ,£n), where ii refers to the index of the fused line segment in slice i, which is 
also the index of the leftmost line segment in the fusing triplet in slice i — I. The set of 
all interaction histories is denoted by (?„. Clearly Gq = and Gn = I mi x . . . x J^^. 
With these preparations the n-th term of the expansion reads 



no+2n „ n n 

TT a{koj,(Joj) / (is(5(t- V'sj)]Texp[-itj(s)fii_i;4(/i;, (Tj,^J] 



as a polynomial of order uq + 2n in the initial fields. Here we have introduced the follow- 
ing shorthands, 

U{s) = ^Sj, i = l,...,n, (2.25) 

i=o 

and 

Vti-jik, a) = -auj{kij + h^j+i + kij+2) - ^{hj) + auj{kij+i) + uj{kij+2) ■ (2.26) 

A„^ contains the 5-functions restricting the integral over k and the sum over a to the 
graph defined by the interaction history i. As can be seen from (12.221) . for non-fusing 
line segments the momentum and the parity are transported unchanged, while the triplet 
of neighboring line segments with parities — , +, + fuses into a line segment with parity 
+ and the triplet of neighboring line segments with parities — , — , + fuses into a line 
segment with parity — . In other words, the middle line conserves parity, while the left 
neighbor carries parity — and the right neighbor parity +. Thus, if the parities in slice n 
are given, then they are already determined for the remaining diagram. At a fusion the 



total momentum is conserved (Kirchhoff s rule). Explicitly, 

2 
j'=0 

X n [^ihj - ki-ij+2)lKj = cr^-l,J+2)]] , (2.27) 



with 1 the symbol for the indicator function. 

The Duhamel expansion of the product (12.231) cut at order A^ is given by 

no N-l 

Y\ a{kj, aj, t) = ^ Tn{t, ki,..., kno, (Ti, ■ ■ ■ , (^no)[^] 
j=l n=0 

+ / dsJ^N{t -s,ki,...,kno,(Ti,..., (Tno) [a{s)] . (2.28) 

As before, in (12.281) a denotes the time field, while a(s) is the time s field as it appears in 
(12.221) . To complete the Feynman diagrams the expansion in (12.281) has still to be averaged 
over some initial state. Note (i) the error term still contains the full time evolution and (ii) 
the operator order for the product appearing in JF„ depends on the particular interaction 
history. 

Given an initial quasifree state of the form (12.171) with bounded W, one can estimate 
roughly the magnitude of a Feynman diagram with respect to this initial state. As an ex- 
ample, let us consider {a{k, t)*a{k', t)) for which uq = 2. At order n of the expansion, 
from the momentum and time integrations one obtains the bound A"t'^c"/ra! with a suit- 
able constant c. The initial state has n\ terms of the same size. Thus a Feynman diagram 
is bounded by A"t"c". However the number of collision histories is also n\, thus yielding 
a zero radius of convergence for the infinite series. 

One possibility to improve the situation would be to extract some extra decay in time. 
A point in case is a Fermi fluid where particles interact only if they are inside a prescribed 
bounded region. Outside this region they do not interact which generates a controllable 
decay. As proved in Il22ll23ll the Duhamel expansion converges in operator norm provided 
the coupling strength satisfies |A| < Aq with small Aq. 

For the kinetic limit no such improved time decay seems to be available and one is left 
with the following tentative program. 

1.) One studies the convergence of each Feynman diagram for rescaled time A^^t, t > 0, 
in the limit A ^ 0. 

2.) Following the pioneering work of Erdos and Yau Il24l in case of the linear Schrodinger 
equation with a weak random potential, one cuts the series at some A-dependent N, the 
generic choice being AA^! = 1, or perhaps A'^A^! = 1 with a suitable choice of k > 0, and 
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tries to control the error term by some other means. In combination with 1.) this would 
allow one to determine the two-point function 

{a{k,X-Hyd{k',X-H)) asA^O. (2.29) 

We will argue that step 1 can be carried out under specified conditions on the dis- 
persion relation u. For step 2, at present our only other help is stationarity. This means 
that we take as initial state (■) the KMS state for H with a suitable choice of the chemical 
potential. Of course, (12.291 ) is then time-independent. Quantities of interest would be two- 
time correlations in equilibrium as {a{k, X^'^t)* a{k')) and, considerably more difficult to 
handle, the number density time correlation {a{k, X~H)*a{k, X~H)a{k')* a{k')) . 

3 The quantum BBGKY hierarchy 

We return to position space. The n-ih. reduced density matrix, at time t is given by the 
expectation 

n 1 

Pn{xi,yi,...,Xn,yn,t) = ( ( JJ a(?/i, t)*) ( JJ a(xj , t)) ) . (3.1) 

Here p„(t) should be regarded as the kernel of a positive operator acting on £2(1^'^)'^"' via 
summation over y. 

For both fermions and bosons Pnit) is symmetric in the arguments (xj,yj), j = 
1, . . . ,n. In addition, for fermions Pn(t) is separately antisymmetric in the arguments 
(xi, . . . , Xn) and (yi, . . . , yn) , while it is symmetric for bosons. To differentiate p„(t) in 
time we use (12.121) . The resulting operator ordering is no longer normal, in general, and 
one has to normal order so as to get a closed evolution equation for the p„(t)'s. Let us 
define, as operators on £2(^'^)^". 

with 

n 

-ffo"Vn(a;i, . . . , a;„) = ^ ^ a{xj - u;)^„(xi, . . . , u;, . . . , x„) (3.3) 

i=i w&Z'i 

and 

n 
^^"Vn(a;i, ...,Xn) = \ ^ V{Xi - Xj)i)n{Xi, • • • , x„) . (3.4) 

We also define 

n 

(C„,„+ip„+i)(xi,...,?/n) = -i XI X] {yi.X3-w)-V{yj-w))pn+i{xu. ..,yn,w,w) . 

(3.5) 
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Then the reduced density matrices satisfy 

^p„(t) = -i[if ("), p„(t)] + \Cn,n+lPn+l{t) , (3.6) 

which is the quantum BBGKY hierarchy. 

On purely mathematical grounds, the unitary evolution extends naturally to the full 
Fock space ^. A state on ^ is given through some positive density matrix S with trgS* = 1. 
We require gauge invariance as [S, e*'^^] = for all 'd E [0, 2n] with N the number 
operator on ^. Thus, denoting by P„ the projection onto the n-particle subspace 3^„, it 
holds 

Sn is a positive operator on ^n and tr^^^^n < 1. Since particles are taken to be indistin- 
guishable, we require Sn to be permutation invariant as an operator on 5„. Let us denote 
by tr[m^„], m < n, the partial trace on ^„ over the tensor product factors with labels m to 
n, in particular trji „] = tr^^. In this general context the n-th reduced density matrix is 
defined through 



oo 



[n + m)\ 



Pn — / i i'^\n+l,n+m]Sn+m (3.8) 

^-^ ml 

as a positive operator on ^„, with the convention tr[„_|_i „]S'„ = Sn- The normalization is 
such that 

tr[i,„]p„ = tT^[SN{N - 1) . . . (iV - n + 1)] , (3.9) 

which is assumed to be bounded by c". Then (|3.8I) can be inverted as 

Sn = / ^ j rtr[„_|_i „_|_miPn+m • (3.10) 

^-^ n\m\ 

m=0 

If tr5[P6)S'] = 1, 9 = ±1, then, using the realization of a{x), a{x)* on Fock space, it is 
easy to check that (13.81) agrees with the definition (13.11) . 
The density matrix S evolves in time through 

Pn,S{t)Pn = 5„„e-^^'"'*^„e'^*"'* . (3.11) 

We insert in the definition (13.81 ) and differentiate with respect to t. The permutation sym- 
metry of the state is preserved in time. This leads to 

-Pnit) = -i[i7("),p„(t)] - iXY,^Iyn+l,n+l\([Vj,n+l,Pn+l{t)]) , (3.12) 

where T^>+i is multiplication by V{xj — Xn+i) as an operator on ^n+i- Clearly, (I3.12h is 



identical to (13.61) . 

The time-integrated version of (13.61) reads 



Pn{t) = e-'^'"'V„e^^'"^* + A f rf.e-'^'"'(*-)C„,.+,p„+,(.)e^^'"'(*-^) . (3.13) 

Jo 
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Of interest is pi(t). Its perturbation series is generated by iterating (13.131) . This will not 
be an expansion in A, since H^''^^ depends itself on A. Expanding exp[—it{HQ + XV^"'^)] 
with respect to A and inserting in (|3.13l) yields a perturbation expansion in A. For bosons 
and fermions it has to be in one-to-one correspondence with the Feynman diagrams of 
Section 2. By construction, for the BBGKY hierarchy the expectation with respect to the 
initial state is over a normal ordered product of operators, while for Feynman diagrams 
normal order does not hold, in general, compare with (12.241) . On the other hand, reading 
Feynman diagrams backwards in time, single line segments can branch only into three 
line segments, while the perturbation expansion of the BBGKY hierarchy in the form 
(|3.13l) consists of n particles interacting amongst themselves and one extra particle added 
through the "collision" Cn,n+i- 

The BBGKY hierarchy is used by Benedetto et al. [[251 |26l HH in their study of the 
kinetic limit. One observes that for the free evolution generated by H^ the kinetic limit, 
space ~ A^^, time ~ A^^, is equivalent to the semiclassical limit. This can be exploited 
by transforming p„ in each of its variables (xj,yj) to a Wigner function. Thereby the 
BBGKY hierarchy turns into a hierarchy of multi-point Wigner functions. The free part 
corresponds to classical particles with kinetic energy uj(k) and the difficulty resides in 
handling the nonlocal "collisions". Benedetto et al. work in the continuum, M'^ instead of 
Z^, and use the quadratic dispersion law uj{k) = k"^. The initial reduced density matrices 
are assumed to be of the factorized form 

n 
PniXi,yi,...,Xn,yn) = JJ Pi (Xj , |/j) , (3.14) 

j=l 

at least asymptotically for small A. (13.141) does not have the antisymmetry required for 

fermions. For bosons the equality P2 = Pi ® Pi forces pi to be a pure state as can be seen 

from 

Remark 3. 1 (Factorization). Using the spectral representation of pi with eigenfunctions 

(pj and eigenvalues \j > 0, symmetry in xi, x^ implies 

oo oo 

^ AiAj0i(a;i)0i(2/i)*0j(x2)0j(2/2)* = ^ AiAj0j(a;i)0i(yi)*0i(a;2)0j(2/2)* • (3.15) 

i,j=l ij=l 

Taking the inner product with 0a(x2) from left and right yields Aq,(1 — Aq.) = 0. <) 

pi defines the scaled Wigner function through 

Wx{r,v,t)= /"t/r/e'^>i(A-V + ir7,A-V-|r/,t). (3.16) 



It is assumed that, for t = 0, 



\miWx{r,v) = WQ{r,v) (3.17) 

A — ^0 
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with Wo{r, v) sufficiently smooth and of rapid decay in both arguments. Benedetto et al. 
prove that the perturbation series for Wx{r,v, A~^t) converges to a limit term by term. 
We describe their limit in Section [51 The convergence imposed in (|3.17l) cannot hold for 
a sequence of Wigner functions coming from a pure state. The factorization (13.141) is 
satisfied only for states which have some support in (1 — P+ — P_)5^ and thus rules out 
bosons and fermions. (13.141 ) is a property characteristic for boltzmannions. 

For bosons and fermions one can switch freely between the BBGKY hierarchy and the 
Duhamel expansion of Section [2l Once we assume the factorization (13.141) . the BBGKY 
hierarchy refers to a larger class of states, not restricted to (P+ + P_) J, and the mapping to 
the Duhamel expansion of Section[2]is lost. Thus in the work of Benedetto et al. some of 
their oscillatory integrals reappear in the Feynman diagrams of the Duhamel expansion. 
But there are still other diagrams. Conversely Benedetto et al. have to consider oscillatory 
integrals which do not correspond to any of the Feynman diagrams studied here. 

4 A comparison with the weakly nonhnear 
Schrodinger equation 

For the Hamiltonian (12.31) we regard a{x) as a complex-valued commutative field and, to 
distinguish, denote it by ^/^ : Z"' ^ C. The classical Hamiltonian functional reads 

H= J2 a{x-y)ij{xrij{y) + lX J^ V{x - y)\ij{x)\'\ij{y)\\ (4.1) 

We set 

H^) = -i^iq, + ip.) (4.2) 

and regard q,j., p^ as canonically conjugate variables. Then the field ^ evolves in time as 

d 
di 



-ij{x,t) = -i 5^ a{x - yMy,t) - iA J^ V{x - y)\ij{y,t)\^ij{x,t) , (4.3) 



which has the same form as (|2.12l) . In particular, using the interaction representation as in 
(12.191 ), the Duhamel expansion for products of the form YYjli i^{^j^ ^j^ ^) is identical to 
(12.281 ) derived for the quantum evolution. The only difference resides in the average over 
the initial state. For the quantum case one has 

no+2n 

(H «(^i'^^)) (4.4) 

with an operator ordering inherited from the Feynman diagram under consideration, while 
in the classical case one has to substitute (14.41) by 

no+2n 

( n ^(^j-'^j))' (4.5) 
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where (■) denotes the average over a suitable initial probability measure on the ^-field. 
Since the ^-field is commutative, the ordering is irrelevant. 

Let us pursue this difference in more detail. For a quasifree state it holds, compare 
with Definition 2.1 and Appendix A, 

2n n 

(JJa(A;j-,o-j)) = XI ^i.'^)W{^i.^^{3)^^-^{j))^i.^An+j)^^-n{n+j))) ■ (4.6) 

j=l 7re'P(2n) j=l 

Here '^{2n) is the set of all pairings of 2n elements with labeling such that, in each factor 
of the product, the operator order is the same as on the left hand side and £:(7r) = 1 
for bosons, e{Tx) = ±1 for fermions depending on whether the permutation induced by 
the pairing is even or odd. Since {a*a*) = = {aa), the average vanishes whenever 
J2j=i ^3 7^ 0- The classical analogue of a quasifree state is a Gaussian measure for which 
the only nonvanishing moments are of the form 

n n 

(ll'^ik.mkn+j)) = Y. U^^ik^y^ikn+.u))) (4-v) 

J = l 7rG-P(n) j = l 

with V{n) denoting the set of all permutations of n elements. Note that (14.71) agrees with 
the bosonic version of (14.61) except for operator ordering. 
For bosons it holds 

{a{kya{k')) = 6{k - k')W{k) , {a{k')a{ky) = 6{k - k'){l + W{k)) (4.8) 

and for fermions 

{a{kya{k')) = 6{k - k')W{k) , {a{k')d{ky) = 6{k - k'){l - W{ky . (4.9) 

If ly is smooth, then so is l+9W{k), and for a subleading Feynman diagram the particular 
operator order makes no difference. On the other hand for a leading diagram one has to 
keep track of the order and the limit will differ classically and quantum mechanically, as it 
should be. In fact, the collision operator of the Boltzmann equation is purely cubic for the 
commutative field while it picks up an additional quadratic piece quantum mechanically 
with a relative sign which depends on the statistics of the particles. 

Remark 4.1 (Rayleigh-Jeans catastrophe). If one replaces Z'^ by W'- and Vl* by —A, then 
(|4.3I) turns into the Hartree equation 

i-Jj{x,t) = -A^{x,t) + X f dyV{x-y)\ij{y,t)\^ij{x,t), (4.10) 

and for a ^-potential into the dispersive nonlinear Schrodinger equation (also called Gross- 
Pitaevskii equation) 

i— V'(a;,t) = -A^(x,t) + X\^{x,t)\^^{x,t) . (4.11) 
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At A = the corresponding equilibrium measure is Gaussian, gauge-invariant, and has 
the covariance 

{iP{ky^{k')) = 6{k-k'){(3{e-fi)y\ ^<0. (4.12) 

For dimension d = 3, the covariance (14.121) is ultraviolet divergent. This is the analogue 
of the classical Ray leigh- Jeans catastrophe for the Maxwell field, in which case the co- 
variance is (5(A; — A;') |A;|^^ (} 



5 The spatially homogeneous Boltzmann-Nordheim equa- 
tion 

We consider an initial state which is quasifree, gauge and translation invariant, and thus 
completely characterized by its two-point function 

{d{kya{k')) = 6{k-k')W{k). (5.1) 

By construction W > and for fermions W < 1 in addition. The dynamics preserves 
gauge and translation invariance. Therefore 

{d{k, tyd{k', t)) = 6{k - k')Wx{k, t) . (5.2) 

As argued above, one expects that, for small A, Wx{k, t) will in approximation be gov- 
erned by a nonlinear transport equation. 

Conjecture 5.1 Under suitable assumptions on uj and on the covariance W in 0. 1 P . it 
holds 

\imWxik,\-h) = W{k,t), (5.3) 

W{k, 0) = W{k), and with this initial condition W{t) satisfies the Boltzmann-Nordheim 
equation 

^W{k,t)=C{W{t)){k). (5.4) 

Here the collision operator is given by 



C{W){ki) = n dk2dksdk46{ki + k2 — k^ — k4)6{uji + uj2 — uj^ — LU4) 

X \V{k2 - fcg) + eV{k2 - k4)\'^{WiW2W3W4 - W1W2W3W4) . (5.5) 

Since the expressions tend to become lengthy we use, here and in what follows, the stan- 
dard shorthand Uj = uj{kj), Wj = W{kj), j = 1, 2, 3, 4. 

Remark 5.2 (Further collision operators). Inserting the definition W = 1 + 9W, the 

16 



Boltzmann-Nordheim collision operator becomes 



C{W){ki) = n dk2dksdk46{ki + k2 ~ k^ — k4)5{uJi + UJ2 — io^ — UJ4) 

x\V{k2 - k^) + eV{k2 - ki)\'^{e{W2W^W^ + WiW^Wi 
-WiW2Wi - W1W2W3) + W^Wi - W1W2) . (5.6) 

In case of the nonlinear Schrodinger equation, d{k) is replaced by the commutative field 
^l'{k). Then (■) is a translation and gauge invariant Gaussian measure and W{k) in (15.11) 
defines its covariance. In this case the collision operator reads 

CNLs(}V){ki) = n dk2dk^dki5{ki + k2 - k^ - k^jS^Ui + U2 - UJ3 - u^) (5.7) 



X \V{k2 - ks) + V{k2 - k^)WW2W^Wi + WiW^^W^ - WiW2W^ - WxW2W- 
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and thus differs from (|5.4I) with 9 = 1 only through the quadratic terms. 

If one imposes the initial condition (|3.14l) corresponding to boltzmannions, then the 
collision operator becomes 



Ccl(W^)(^i) = 2n dk2dk^dk46{ki + k2 — k^ — k4)6{uji + UJ2 — ujs — UJ4). 

x\V{k2-ks)\\W3W4-WiW2). (5.8) 

In their set-up Benedetto et al. prove (15.81 ) in the sense that the perturbation series gen- 
erated by the BBGKY hierarchy (13.61 ) converges term by term to the perturbation series 
generated by (15.41 ) with collision operator Cql- Since they do not renormalize the dis- 
persion as in (16.111) below, they have to impose that V{0) = 0. One recognizes (|5.8I) as 
the classical Boltzmann equation with the Bom approximation to the differential cross 
section. Thus on the kinetic level boltzmannions behave like classical point particles. 
Note that for bosons, 9 = 1, 

C = Cnls + Ccl, (5.9) 

to say, adding the collision operators for classical waves and classical particles yields the 
quantum mechanical collision operator. It is surprising that quantizing either the nonlinear 
Schrodinger equation or classical point particles results in such a small modification on 
the level of the kinetic equation. (} 

Remark 5.3 (Spatially inhomogeneous Boltzmann equation). If one adds in (15.41) . (15.51) 
the spatial variation, then the Boltzmann-Nordheim equation becomes 

^W{r, k, t) + ^VkOo{k) ■ VrW{r, k, t) = C{W{r, -, t)) {k) , (5.10) 
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where our notation is supposed to indicate that C acts on the argument k at fixed r, t. Of 
course, the same holds for C replaced by Cnls or Cql- (|5.10I) can be interpreted as com- 
ing from the motion of classical particles with kinetic energy uj(k). Collisions between 
particles are implicitly defined through the conservation of energy and momentum, 

tui + iU2 = ujs + UJ4 , ki + k2 = ks + k4. (5.11) 

The collision rule thus depends on the particular form of u and can be very counterin- 
tuitive when viewed from the perspective of potential scattering of mechanical particles. 
In the case of the special dispersion relation uj{k) = k'^/2 on M^, energy and momentum 
conservation can be parameterized in the form 

ks = ki — u ■ {ki — k2)Cj , k^ = k2 + Cj ■ {ki — k2)Cj , (5.12) 

with \Cj\ = 1, i.e., u E S^. Then the collision operator acquires the more conventional 
form 

C{W)ik,)= [ dk2 [ duj\uJ-{k,-k2)\\V{iU-ik,-k2)iu) 

+9V{ki -k2-uj-{ki- k2)Cj)\^[WiW2W^Wi - WiW2W^Wi) . (5.13) 

For a rotational symmetric potential, V{k) = Vr{\k\), the collision cross section simplifies 
to 

\U ■ {k, - k2)\mUJ ■ (k, - k2)\) + 9Vr{m - hf - {Cj ■ {k, - k2)fY'^) ^ . (5.14) 

One notes that for a smooth potential the decay is exponential and even for hard spheres 
the decay is proportional to \ki — k2\^'^. Thus at high energies there is only little scatter- 
ing. <> 

Remark 5.4 (History). Equations (15.41) . (|5.5I) were first written down by Nordheim [HI, 
where he had in mind the true quantum mechanical scattering cross section, rather than 
only its Born approximation. In 1929 Peierls [2] studied lattice vibrations with small non- 
linearity both classically and quantized. For this particular weakly nonlinear wave equa- 
tion he derives (|5.4I) with the analogue of the collision term (15.71) . Later on it was realized 
that Peierls' ideas apply to a more general class of weakly nonlinear wave equations, e.g. 
see [|281 . For quantized lattice vibrations Peierls uses Fermi's golden rule and arrives at 
(15.41) with the analogue of the collision operator (|5.5I) for 9 = 1. The terminology is not 
uniform. In kinetic theory the name Uehling-Uhlenbeck seems to be most frequent be- 
cause they studied the equation in their pioneering work [3]. For phonon transport Peierls 
or Boltzmann-Peierls is used. For dilute Rose gases Boltzmann-Nordheim seems to be 
rather established. We follow this latter convention for reasons of priority. (} 

To approach the Conjecture one expands {d{k, t)*d{k', t)) in the Duhamel series (|2.28l) 
up to some N depending suitably on A. Currently there seems to be no good idea of how 
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to control the error term. This leaves one with the program of the term by term conver- 
gence, which to some extent will be explained in Section [6l An important improvement 
as regards to Section [21 is to renormalize the bare dispersion uj to 

u:\ki, t) = uj{ki) + \ [ dk2Wx{k2, t) (f (0) + eV{ki - k2)) . (5.15) 

In contrast to Section [61 in the present context to'^ is time-dependent through the itself 
unknown Wx{t). 

By mass conservation 



dkiW{ki,t)= dkiW{ki,0). (5.16) 

jd Jjd 

Thus the term proportional to V{0) is in fact constant. For the second term one has the 
trivial estimate 

A f dk2Wx{k2,t)9V{ki - Aia) <\f dkiW{ki,0) V \V{x)\ . (5.17) 

These observations leave us with two choices: If V{k) = V{0), then the fermions become 
noninteracting and one is left with lattice bosons interacting through a quartic on-site 
potential, cu is renormalized by a constant proportional to A which is easily taken care off. 
If V{k) depends on k, the situation is more complicated. Because of the convolution, u;^ 
depends smoothly on fci, but it could be rapidly oscillating in t. Whether the bound (15.171 ) 
suffices to ensure the £3 dispersivity, the constructive interference, and the crossing bound 
of Section [6lremains to be investigated. 

In this section we discuss a more modest step, namely the leading part of the main 
term. For this purpose we make the following definitions. 

Definition 5.5 (pairing property). Let us consider a Feynman diagram of even order, n 
even, where the integration over all momentum S-functions has been carried out. The 
Feynman diagram satisfies the pairing rule if for every even time slice the momenta are 
paired, i.e., to each line segment with momentum k and parity a there exists, in the same 
even time slice another line segment with momentum —k and parity —a. 

Definition 5.6 (leading diagrams). A Feynman diagram is leading, if it satisfies the pairing 
property and if it does not contain the factor 1^(0). 

Note that in the 0-th time slice [0, sq] the pairing is induced by the initial state. Thus the 
structure of a leading diagram can be obtained by iteration. 

So let us assume that for the time slice with even label i all line segments are paired. 
If there are 2n line segments, they carry momenta cTjkj, j = 1, . . . ,n, aj = ±1. kj 
and —kj are paired. At the end of time slice i three neighboring lines fuse to the line 
segment with momentum a'k'. If two momenta would be paired, the diagram necessarily 
contains the factor V{0). Thus we label neighboring lines by three distinct momenta o-jkj, 
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Figure 2: A leading Feynman diagram with uq = 4, n = 4, and interaction history 
(3,3,5,1). Up (down) arrows correspond to parity +1 (—1). The pairing from the initial 
state is indicated, where right (left) pointing arrow stands for the order {a*a) ({aa*)). 



j = 1,2,3. In the next step there can be two cases which are both displayed in Fig.[2l Case 
a) corresponds to the 3"* and 4* fusion, while case b) corresponds to the P' and 2"'' fusion. 
In case a) the line segment a'k' does not participate in the fusion at the end of time slice 
i + 1. The fusing triplet has momenta aT^(^j)k^(^j), j = 1, 2, 3, and the momentum of the 
fused line segment is a"k". By the same argument as above, the 7r(j)'s must be distinct. 
To have pairing in time slice i + 2 requires k' = k", a' = —a". In turn this is possible 
only if {1, 2, 3} = {7r(l), 7r(2), vr(3)}. In case b) the line segment a'k' participates in the 
fusion at the end of time slice i + 1. The fusing momenta are a'k', o",r(i)^7r(i), cr7r(2)fc7r(2)- 
The fused momentum is a"k". In slice i + 2, a"k" must be paired with, say, cr^(3)A;^(3). 
There will be a factor 1^(0) unless all three 7r(j)'s are distinct. The pairing of a"k" and 
crn{3)K{3) is possible only if {1, 2, 3} = {7r(l), 7r(2), 7r(3)}. 

We conclude that for a leading Feynman diagram there must be three paired line seg- 
ments in slice i which connect through two fusions to a single paired line segment in slice 
i + 2, i = 0,2, ... ,n,n even, compare with Fig. [21 This property allows us to represent a 
leading diagram through a contracted diagram which we explain next, see Fig. [31 

In a contracted diagram we draw only the even time slices and each pair as a single 
line segment. A line segment thus carries a momentum k. But we still have to distinguish 
the relative order within the pair. If in the original Feynman diagram the order from left 
to right is — , +, then the line carries the order parity r = +1, while for the order +, — the 
order parity is r = — 1. In the contracted diagram, at the end of each time slice a triplet 
of neighboring line segments fuses into a single line segment. In our case, if the leading 
Feynman diagram has order 2n, then each contracted diagram has n collisions and there 
are (2n + l)!/2"n! contracted diagrams. 

To compute the vertex strength for the contracted diagram one has to sum over all 
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Figure 3: Contraction of the diagram from Figure 2. The order parity, r, of each con- 
tracted line is indicated. 

Feynman diagrams which fuse 3 pairs at even time slice into the single pair at time slice 
i + 2. There are 16 such diagrams. Their sum yields the vertex strength given below in 

dm. 

To write down the integral corresponding to a contracted diagram, it is useful to intro- 
duce the "correlation" functions p„ : (T'^x {— 1,1})"^C and to define time-dependent 



collision operators Cj,„+2(i) with j 



n. We set 



{Cj.a+2{t)pn+2) {h, n, . . . , K, r„) = 2 / dg2dg^dg^5{kj + 92-93- 94) 

J{Td)3 

xcos{{iu{kj) +uj{92) -io{9s) - uj{94))t)V{92 - 93){Vi92 - 93) +dV{92 -94)) 



X (p„+2(- ..,92,-Tj,..., 93, Tj, 94, Tj) + 6'p„+2(. • • , kj, Tj , 
-pn+2{- ..,kj,Tj,..., 92, 1, 93, -1) - 0Pn+2( 

which acts on the arguments j, n + 1, n + 2 of Pn+2, and 



, kj , Tj , 



Let 



Cn+2{t) — 7 ^Cj,n+2{t) 

M^ik,,t)5iko-k^) 



,5'3,1,5'4,1) 
,^2,1,^4,1)), (5.18) 



(5.19) 
(5.20) 



be the sum over all contracted diagrams at order n with time span t. Then 

M^^ik^,X-H) 

P „ n+l n 

= A'" / dt ds6{y^tj + Vs, - X-h){C3{Si) . ..C2n+l{Sn)p2n+l){kl, 1) 

= dt ds{C3{Si) ...C2n+l{Sn)P2n+l){kl,l) . (5.21) 

In our particular case the input function is given by 

n 

pn{ki, Ti, . . . , fc„, r„) = JJ W{kj, Tj) (5.22) 

3=1 
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with 

W{k,l) = W{k), W{k,-l) = l + eW{k) = W{k). (5.23) 

To prove the limit A — i> one needs a first assumption on to. 
Assumption Al (£3 dispersivity). Let us define 



Then there exist c > 0, 5 > such that 



I rfA;e-'*^('^)e'2^^-^ (5.24) 



Y.\pt{x)\'<c{tr'~' (5.25) 



with the shorthand (t) = vT+t^. 

As in [[T5I one proves the following bound. 

Proposition 5.7 Let \\pn\\i < (co)", where || ■ ||i is the £i(Z'^)®" norm in position space. 



Then there exists a constant c such that 

r"00 



/ rfs sup I (C3(si)...C2„+i(s„)p2n+i)(A;i,l) I < c^'n!. (5.26) 

Note that C2n+i{s) has 2n — \ terms of equal size. The n\ in (15.261) thus results from the 
product of collision operators. 

With this bound one can introduce the collision operator 



{Cp^{ki, t) = 2n dk2dk-i6{uj{ki) + uj{k2) - uj{k-i) - tu{ki + k2 - k^)) 

J{T'i)2 

xV{k2 - k-i) {V{k2 - k^) + eV{k2 - k^)) {^^2, -r, k^, r, k^ + k2 - k^, r) 
+9ps{ki, r, k2, 1, ki + k2- k^, 1) - ^3(^1, r, k2, 1, k^, -1) 
-eMh,r, /C2, 1, ki + k2- ks, 1)) , (5.27) 

where the ^-function is defined through the limit 

5(n) = \im-—^. (5.28) 

As before, from C we construct Cj „+2, which is C now acting on the variables kj, kn+i, 

kn+2, and 

n 
Cn+2 = 2_^ Cj,n+2 ■ (5.29) 

Under the conditions of Proposition 5.7 one concludes that 

\im Mi{k^,\-H) = -e{C^...C2n+lP2n+l){kl^) (5.30) 

A — >0 Tl\ 
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uniformly in ki and 



sup |(C3 . . .C2„+iP2n+i)(fci, 1)1 < d'nl . (5.31) 



Thus the sum over n reads 

^(^' ^) = Zl ^(^3 . . . C2n+lP2n+l){k, 1) . (5.32) 



oo ^, 



n! 

n=0 



By (15.261) there exists then a Iq, such that the sum converges provided 

0<t<to- (5.33) 

If of interest, the time to can be computed explicitly. Up to numerical factors of order 1 , 
to is proportional to 

dtY,\Pt{^)?- (5.34) 

Because of the particular initial conditions (15.221) . (15.231) . the limit in (|5.32l) can be 
written more concisely. One recognizes W{k, t) as the power series solution of the non- 
linear Boltzmann-Nordheim equation (15.41) with collision operator 

C{W){ki) = 27r / dk2dk'idk/^5{ki + k2 — k^ — k4)6{uji + uj2 — ^3 — ^4) 

J{J'i)3 

xvik2 - ks){vik2 - k3) + evik2 - k^)) 

X {W2WsW4 + eWiWsW^ - W1W2W3 - eWiW2W^) , (5.35) 

which is identical to (15.51) . 

6 Two-point time correlations in thermal equilibrium 

6.1 Set-up. In this and the following section we consider a quantum fluid in thermal 
equilibrium. We fix some inverse temperature /3 > and a chemical potential fi eM. such 
that 

u{k)-ii>0 (6.1) 

for all k E T'^. (■)/3,a denotes the expectation value with respect to the /5-KMS state for 
the hamiltonian H — fiN, N the number operator. The two-point function is denoted by 

{d{kyd{k'))p,x = 6{k - k')wl^{k) , 

(a(A;)a(A;')*)ftA = S{k - k')Wl^{k) = S{k - k') (l + 9Wl^{k)) . (6.2) 

In the limit A ^ 0, the state (■)/3,o is quasifree with 

{d{k)*d{k'))p^o = 6{k - k')W^{k) . (6.3) 
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Since Wj^ is smooth, the two-point function decays exponentially and by quasifreeness 
all fully truncated correlation functions of order greater than two vanish. 

In this section we study the effective propagator as the most basic two-point function, 
namely 

{a{kya{k', t))f,,x = {a{k)*a{k'))p,xCxik, t) , (6.4) 

which defines C\{k, t) since both sides of (|6.4I) are proportional to 5{k — k'). More com- 
plicated equilibrium time correlation functions will be discussed in Section 7. Because of 
the a-term, clearly, C\ will contain the oscillatory factor e~"^(^^* varying on time scale 1 . 
To first order in the Duhamel expansion one has 

e'-('^^)*(a(A;o)*a(A;i,t))^,;, = (a(A;o)*a(A;i))/3,A 

— iA I ds I dk2dksdk45{ki + k2 — k^ ~ k4)V{k2 — k^) 

Jo J{Ti)3 

X exp[is(co'i + UJ2- UJ3- UJ4)] {d{ko)* d{k2)* d{ks)d{k4)) fs^x + C(A^) . (6.5) 
We write 

{d{ko)*d{k2)*d{k3)d{ki))i3^x = {d{ko)*d{k^))f3^x{d{k2)*d{k3))0^x 
+0{d{ko)*d{k3))f^^x{d{k2y*d{k4))i3^x + {a,{kQ)*a{k2)*d{k^)a{ki))1y^ . (6.6) 

The truncated part is 0{\) and, when first integrated over momenta as in (16.51) . is abso- 
lutely integrable in s. Therefore this contribution will be part of the error term. The first 
two terms inserted in (|6.5I) yield —iXt{d{ko)*d{ki))i3^xRxiki), where 

Rx{k,) = [ dk2W'x{k2) {V{0) + eV{ki - k2)) . (6.7) 

Our computation suggests that Cx has a second oscillatory factor of the form 

varying on the time scale A^^ The kinetic time scale is order A^^. Thus it suffices to 
expand up to A as 

Rx{k^) = Ro{k^) + Ai?i(A;i) + ^(A^) (6.9) 



with 



Ro{k,)= [ dk2Wl^{k2){V{0) + eV{ki-k2)), 

Riih) = -P idk2Wl{k2)wl{k2) (t/(o) + ev{k, - k2)) 
dk^wf,{k^){v{S)) + ev{k^ - ks)) . 



Ijd 

X 



(6.10) 
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Potentially the oscillatory term (16.81) could be dangerous because in an expansion in 
A it may mask the kinetic terms. As we will see, fortunately, Rx can be absorbed by 
renormalizing cu to 

uj^ = uj + XRx. (6.11) 

The second order term of the Duhamel expansion will not be written out explicitly. It 
consists of 18 diagrams, 12 of which combine to 

- iAV/?o(A:)' . (6.12) 

The remaining 6 diagrams sum up to 

- 6{k - k')W^^{k)u{k)XH (6.13) 

valid for large t. The decay coefficient u is obtained as 

z/(A;i) = — / dt dk2dk^dk/^5{ki + k2 — k^ — k^) exp[it(c(Ji + uj2 — to^ — cu^)] 

Jo JiTdyi 

y<v{k2 - ks){v{k2 - k,) + ev{k2 - k,)){wl,wl, - wl^w^, - Qw^^wl^) 

(6.14) 
with the shorthand W^ ■ = W^f^{kj). For the real part of v one obtains 

'^v{k{) = IT dk2dkj,dki5{ki + k2 — k^ — k4)S{uJi + uj2 — co-^ — cu^) 

xl{V{k2 - ks) + eV{k2 - k^))\wl,)-'Wl2Wl^Wl^. (6.15) 

In particular 

^u{ki)>0. (6.16) 

On the basis of this second order expansion the obvious conjecture is that 

Cx{k,t) = exp [ - {iuj^ + X'^u{k))t] (6.17) 

for t > and t = C(A^^), the case t < following from time reversal as 

Cx{k,ty = Cx{k,-t). (6.18) 

Our goal is to prove (16.171 ). For this purpose we have to assume cluster properties of 
the equilibrium state and the decay of certain oscillatory integrals. 

6.2 li-clustering of fully truncated correlation functions. For the KMS state (■)/3,a we 
consider the fully truncated correlation functions denoted by (\Xj=i^{^3i^3)%,\- We 
refer to Appendix A for their definition. They vanish whenever YTj=i ^j 7^ ^^'^^ ^^ 
proved in Appendix A, for n> A they do not depend on the operator ordering except for 
an overall sign. 
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Assumption A2 (£i-clustering). Let /5 > and ^ satisfy rt<5JD . There exists Aq > and 
Co > independent ofn such that for < A < Aq. and alln > 4 one has the bound 

n 

^ S.J{l[a{x„aj))l^ <A(corn!. (6.19) 






In addition, 



J2 I («(0)*a(x));3,A - (a(0)*a(x));3| < \2{coY ■ (6-20) 



Ginibre [jTTll studies £i-clustering for low density quantum gases in the continuum, 
i.e., for position space M^ and dispersion relation uj(k) = k"^. It is rather likely that 
his analysis could be carried through also for lattice gases. The L^ -bound by Ginibre is 
based on an expansion with respect to the fugacity, hence the small coupling regime is not 
optimally covered. In particular, the prefactor A in (16.191 ) cannot be deduced from [fTTl . It 
seems that extra work is needed. 

6.3 Oscillatory integrals. The ^3-dispersivity has been stated already in Assumption Al. 
In addition we need an assumption which controls the constructive interference between 
two frequencies. 

Assumption A3 (constructive interference). There exists a set M^™s q f c( (^Qfisisting 
of a union of a finite number of closed, one -dimensional, smooth submanifolds, and a 
constant C such that for all t G M, A;o ^ T'^, and a G {±1}, 






^ C{ty' r6 2n 



where d{ko, M'™^) is the distance of ko from M''™^. 

Remark 6.1 (Dimension). Assumption A3 allows us to cut out a small tube around each 
curve. If one would cut out too much, e.g. two-dimensional surfaces, this will show up 
in other parts of the proof. For this reason Assumption A3 as stated requires in addition 
d > 4. If (|6.21l) would hold for M*^™^ merely a collection of a finite number of points, 
then we could accommodate d > 3. (} 

Next we need a mechanism which allows to distinguish between leading and sublead- 
ing diagrams. For the linear Schrodinger equation with a random potential, this mecha- 
nism has been identified in [24] . In the graphical representation developed in this article, 
it corresponds to the crossing of two edges. Our condition is the natural generalization of 
the crossing estimate in Il24l . 

Assumption A4 (crossing bounds). Define for to, ^1,^2 G M, mi, U2 G T'^, and x G Z'^, 

K{x;to,ti,t2,uuU2) = [ cifce'2""-'=e-'(*°"('^)+*^"('=+"^)+*^"('=+"^)) . (6.22) 

We assume that there is a measurable function F^"^ : T^ x ]R_|_ — > [0, 00] so that constants 
< 7 < 1, Ci, C2for the following bounds can be found. 
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(i) For any Ui G T*^, cXi G {±1}, i = 1, 2, 3, and < /3 < 1, the following bounds are 
satisfied: 



||2 



/ dsids2e ^^''''+''''^||i^(si + S2,criS2, Class, Ml, M2)||3|bsi- 

<(3'<-^F'\u2-ur,(3), (6.23) 

/ rfsirfs2 e-^(l^^l+l^^l) n \\K{si + 52, a,S2, 0, n„ 0)||3 

< P^-^F^'{un; P), for any n G {1, 2, 3} . (6.24) 

(ii) For all < ( < (q we have 

[ dkF^'{k-C) <Cl(lnC)'^ (6.25) 

and if also u,ko G T'^, a E R, a E {±1}, and n G {1,2,3}, and we denote k = 
{ki, /c2, ko — ki — ^2), then 

[ dk,dk2F^^{k^ + u;Ch ^.} ,^.^1 < Cl(lnC)^+^^ (6.26) 

where Q : (T'^)^ x {±1} -^ M is defined by 

n{k, a) = uj{k^) - uo{ki) + a{uj{k2) - uj{ki + k2 + k^)) . (6.27) 

Remark 6.2 (Reduction). The reader may wonder whether the stated assumptions can 
be proved for a specified class of cij's. Let us first emphasize that the reduction of the 
many, very high-dimensional oscillatory integrals to a few low dimensional oscillatory 
integrals involving only uj is already a big step in the right direction. It is also clear that 
the reduction cannot be pushed any further. The £3 dispersivity is needed already to define 
the limit equation. Constructive interference occurs always at fco = 0. Thus some control 
of the phenomenon is required. The crossing bound reflects the mechanism how only a 
few diagrams survive in the limit. At present, A4 is one version which works, but its 
optimal form could eventually look differently. In any case. Assumptions Al, A3, and 
A4 rely on results from a disjoint mathematical discipline. In this sense, the situation is 
similar to A2. Clustering is proved within rigorous statistical mechanics. In our context it 
is a necessary input stated in a form which is regarded as obvious by the experts. Al can 
be proved by stationary phase methods and holds generically for dimension rf > 3. We 
are not aware that A3 has ever been studied systematically, but hope that our work might 
serve as a motivation. A4 in the somewhat simpler context of the random Schrodinger 
equation has been proved for nearest neighbor couplings in [|29ll and is investigated more 
systematically in (SOI [HI • 

Remark 6.3 (Example). For on-site and nearest neighbor couplings only, i.e., a{x) = 
for \x\ > 1, a{x) = «! for |x| = 1, and arbitrary «(0), Al is proved in [,19ll for d > 3 and 
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A3, A4 are proved in [f TSl for d > A. The choice of the constants is 7 = |, C2 = 0, and 
the function F^^ is taken to be 



d 



F-{u-X) = Cll . J (6.28) 

with a certain constant C depending only on d and to. (} 

6.4 Main result. To formulate our main result it is convenient to integrate the field against 
a test function f E £2- Interpreting (-, ■) as inner product, we set 

(/, a)=J2 f{^r<^) , (/, a)* = ^ f{x)a{xr (6.29) 

and correspondingly for our other conventions as (/, a(t)) and (/, a(cr, t)). In particular 



(/,a)= / dkf{kyd{k) = {f,a). (6.30) 

Note that by Schwarz inequality for operators and stationarity 

|((/i,a)*(/2,a(t))>^^j' < ((/i,a)*(A,a)>^^,((/2,a)*(/2,a)>^^, 



dkWl^{k)\h{k)\' dkW;^^{k)\Uk)\- (6.31) 

Hence the quadratic form (/i,/2) 1— > ((/i, a)*(/2, a(t)))g^ is uniformly bounded on 
£2 X ^2- Our main result concerns the behavior of this quadratic form on the kinetic time 
scale. 

Theorem 6.4 Let {■)i3,x satisfy Assumption A2, let uo satisfy Assumptions Al, A3, andA4, 
and let d> A. Then there exists to > such that for <t < to one has 

lim ((/2,a)*(exp[-itA-V]/i,a(A-2t))> = / dkWl,{k)j\{kye~'''-'^'f2{k) 

(6.32) 

/ora///i,/2G£2W. 

Remark 6.5 (Short kinetic time). The restriction to a short kinetic time, t < to, has no 
physical significance. The decay is expected to hold for arbitrary to and most likely for 
even longer times. The smallness of to merely reflects that in our context it is difficult 
to properly bound the error terms. Technically it arises because in the error term one 
obtains high order diagrams which are close to the spatially homogeneous case as studied 
in Section [5l One needs that the leading part of the main term is small. Since the leading 
diagrams are bounded as (t/to)^, one has to require t < to. <) 

6.5 Link to the nonlinear Schrodinger equation as discussed in / f75l/ . [[TSlI is written in 
such a way that the operator ordering, as of relevance in our context, is already respected. 
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Thus, with the proper reinterpretation, most formulas hold also for the quantum evolution. 
In particular, the analysis of the Feynman diagrams carries over verbatim. The Duhamel 
expansion (12.281) is preliminary and one still needs three modifications to reach the starting 
point of [15]. 

— partial time integration. This refers to the time integrations and thus remains valid in 
the quantum mechanical context. 

— insertion of the cutoff functions <I>q,$^. This refers to the fc-integrations and thus 
remains valid in the quantum mechanical context. 

— removal of fast oscillations. As argued above the dispersion relation uj is renormalized 
to u;'*'. Thus, as a modification of (12.191 ). we define 

a(fc, 1, t) = e''"'('=)*a(A:, t) , a{k, -1, t) = e"'^^ ^^^' a{-k , t)* . (6.33) 

In addition we introduce the pair truncation 

V{a{ki, -l)a{k2, a)a{ks, 1)) 

= a{ki, -l)a(/c2, cr)a{k3, 1) - {{a{ki, -l)a{k2, (r))p^xa{k3, 1) 
+9{a{ki, -l)a(/c3, 1))i3,\a{k2, a) + {a{k2, <7)d{k3, l))f3^xa{ki, -1)) , (6.34) 

a = ±1 and omitting the t argument. Then a{ki, a, t) satisfies the evolution equation 

—a{ki, a, t) = —iXa / dk2dk'idki5{ki — ^2 — fca — ^4) 

dt J(T<*)3 

xi((l + a)V{k2 + k;) + (1 - a)V{k3 + k^)) 
^g-itn(fc,<x) i^A^^^^ A;3, ki)d{k2, -1, t)a(A;3, a, t)d{ki, 1, t) 
+^^ik2, k^, k^)P[d{k2, -1, t)a(/c3, a, t)d{k^, 1, t)] } , (6.35) 

where 

n{k, a) = -auj^iki) - Lj^{k2) + au^iks) + uj^{ki) . (6.36) 

(|6.35l) agrees with the corresponding formula of [flSl with the only modification consisting 
of a general interaction potential. In [15J we consider the case V{k) = V{0). Therefore 
the factor in (16.351) containing V is replaced by a constant and uj'^{k) = uj{k) + XRq, 
where Rq depends on A but not on k. In frequency differences Rq thus drops out. 

We are confident that the estimates of the oscillatory integrals proved in ifTSl still hold 
for the case under consideration here. Of course, a simple estimate as 



e 



-iu{k)t _^-iu\k)t^^^m^ (6.37) 



uniformly in k is too crude. Rather each oscillating integral has to be estimated with an 
uj^ depending weakly on A. This step remains to be carried out. 

The classical average E in [15^ becomes the quantum average (■)/3,a over the KMS 
state. With the assumption A2, higher cumulants do not contribute. For the quasifree 
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part one has either W or W as time t = input. Since both functions are smooth, the 
estimates from the commutative case remain valid. 

6.6 Error terms. We adopt the notation from IfTSl . There are three error terms of the same 
structure, namely 

((/i,a(0))* / ds{h.X^{t,s)[a{s)])) (6.38) 

Jo 

with X any of Q, Z, or A. We use the Schwarz inequality for operators A, B according 

to which 

l(^*i?WP<(|AP)/3,A(]i?H,,A (6.39) 

with the shorthand A*A= \A\^. Then 



t3,X 




< ( / ds\{{h,d{o)){kMs,t)[d{s)])) 



I3,X\ 



<i(|(/i,a(0))|\,y^ ds{\{h,X^{t,s)[d{s)])\')^^^. (6.40) 

Xn{t, s) is a monomial of order n in the factors a{k, a, s) which differ by a phase factor 
from a{k, s), a(k, s)*. The latter field operators are invariant for (■)/3,a and, at the expense 
of a phase factor, d{s) in (16.401) can be replaced by a(0). Thus, at given cutoff N, the 
two-point function of (16.321) is reduced to a main term and an error term, which both 
involve only the free time evolution and are monomials of the time fields a averaged 
with respect to the /3-KMS state (■)/3,a- 

6. 7 Convergence of the leading part of the main term. We refer to Section[5]for the general 
structure. For the case under consideration the Feynman diagrams are constrained. In 
(16.321) the oscillating term is combined with (/2, a)*. The Feynman diagram has one line 
segment [0,t] with parity —1, which is constrained not to fuse at all. In the time slice 
[t — s„,t] there is only one further line segment. Thus ng = 2. A leading diagram is 
still defined by the pairing property for even time slices and the absence of the factor 
y(0). The line corresponding to (/2, a)* has momentum ki. The collision operator acting 
on k2, ■ ■ ■ , kn+2 is defined as in (15.181) . Because of the constraint the collision operator 
acting on ki reads 

{<^i,n+2it)Pn+2) (^1, Ti, . . . , /c„, r„) = 2 / dg2dg^dgi5{ki +92-93- 9i) 

^git(<^(fci)+<^(32)-<^(93)-<^(94))y(^2 _ k3){V{k2 - ks) + eV{k2 - k^)) 

X {9pn+2iki, Ti, . . . , 93, 1, 94, 1) - pn+liki, Ti, . . . , 92, 1, 5-3, "1) 

-9pn+2iki, Ti, . . . , 5f2, 1, 5-4, 1)) ■ (6.41) 

The initial conditions are 

n 
Pn{ki, Ti,...,kn,Tn) = Y\_ ^f^i^j, Tj) . (6.42) 

i=i 
30 



We set 

n 
C-+2it) = C{,n+2it) + Y. '^^•'«+2(t) . (6.43) 

i=2 

Let us define A^^(t) as the sum of all leading Feynman diagrams at order n. With the 
above notation it is given by 

„ „ n+l n „ 

>JnW = A'"/ dt dsdiTtj + Ts — t) dhAihrUh) 

J(M+)"+i J(R+)" j^^ j^^ J-Fd 

{C^is,) . . .C2-+i(s„)p2„+i)(A;i, 1) . (6.44) 

Proposition 5.7 remains valid. Thus one can pass to the limit as 

limA^^(A-^t) = ^r / dhf,ihyf2{h){C^ . . .C2;+iP2n+i)(A:i, 1) • (6.45) 
A^o n! Jfd 

Let ps be given by ( 16.421 ). Then 

(Ci,3(s)p3) {ki, n) = — / dk2dk3dh6{ki + k2 - k^ - k^) 

x{siYiVLs)^-\e^^ - l){V{k2 - ks) + e{V{k2 - k^))^ 

xWfsik,, n)Wp{k2, n)Wf,{k,, -n)Wf,{k,, -n) , (6.46) 

with n = LUi + UJ2 — co-i — UJ4. Therefore 

/ ds{C^,s{s)p3)ih,r^)=0. (6.47) 

Returning to (16.451) we conclude that when acting with Cg^+i, where j = 1, . . . ,n, only 
its first summand, i.e., C^2j+i^ contributes. Hence, there exists to such that for < t < to 
one has 

lim VA<^(A-^t) = Y^-(-lY / dkhiky f2{k)Wl{k)v{kr . (6.48) 

n=0 n=0 -^ ^ 

Remark 6.6 (Radius of convergence). The integral in (|6.46l) decays at least as s^^ . There- 
fore the error in (|6.47l) is order A^*^. The leading diagrams at order n have a single term 
which yields the non-zero limit in (16.481) . while the remaining (n! — 1) terms are of order 
A^'^. Therefore one can cut the series at order N with A^'^A^! = 1. Since the series of con- 
tributing terms has an infinite radius of convergence, the restriction to bounded to can be 
lifted and (|6.48l) holds in fact for all t. On the other hand the error term can be controlled 
only for \t\ < to and the extra effort will not improve Theorem 6.3. <0> 
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7 Time correlations for the number density 

Physically of great interest are the density current and energy current time correlations. 
They differ in two respects from the correlation function studied in the previous section: 
they are four-point functions and, more importantly, involve a spatial summation. Slightly 
formal, these correlation functions are particular cases of the number density time corre- 
lations in equilibrium, i.e., of {a{k, t)*a{k, t)a{k', 0)*a{k', 0))i3,x. 

Let 1] : Z^ -^ C such that t](x)* = r]{—x) and r] has bounded support. We define the 
energy-like, resp. number-like, observable 

H^= Yl vi^-yHxYaiy) (7.1) 



and represent it by a sum of local terms as 



W 5 



^^ = lY (a(w;)*?7(w - y)a{y) + a{y)*r]{y - w)a{w)) . (7.2) 

For H"^ = i^har one has r)(A;) = uj{k), for the energy current f]{k) = {SJ k^{k))uj{k), and, 
correspondingly, T){k) = 1, r](k) = Vk^{k) for number and number current. 

As before (•)/3,a is the /9-KMS state for H — fiN at infinite volume with /x satisfying 
(16.11) and we set (■)^ o = (■)/3- It is also convenient to introduce the Kubo inner product 
defined through 

{{A, 5))^,, = p-' f dp'{{A*e^^'^Be^'^)f,,, - (A*)^,,(S)^,,) . (7.3) 

Jo 

The fluctuations of H^ in a large box A are given by 

weA 
and the quantity of interest is the time-displaced covariance of (|7.4I) . which reads 

Jim ((a, am,,x = J2 ((^- ^0 (^)))aa = cut) . a. 5) 

toGZ'* 

In particular, at t = 0, 

limlim((a,a))AA= E«^-^o))/5= / dk\f,{k)\'W',{k)W;{k) . (7.6) 

Note that in the limit A -^ the Kubo and standard inner product coincide for the observ- 
ables under consideration. 
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As before the kinetic limit provides information on C^(A ^t) for small A. There is 
a simple formal argument how to guess this correlation by building on the results from 
Sections Let {■)f3,xie) be the (3-¥MS state for H - fiN + {e/l3)H'^. Then 

We now perform first the limit A ^ and then the limit £ ^ 0. For time t = one has 

lim(a(A;)*a(A;'))/3,A(£) = 5{k - k')W'{k) , (7.8) 

A — >\j 

where 

W'{k) = (e[/3(^W-M)+£'?W] _ gyi ^ (7_9) 

which is well-defined for e sufficiently small. Now, by Section [51 

lim(a(A;, X-Hya{k', X'H))f3^x{e) = S{k - k')W%k, t) , (7.10) 

A — >{) 

and W^i^k, t) solves (|54l) with initial condition W'^{k). It follows that 

\im{H^,{X-h))p,x{e)= [ dkfj{k)W%k,t) . (7.11) 

We linearize the collision operator as 

C{W'p+ef) = EAf + Oie'). (7.12) 

It is convenient to introduce the multiplication operator U/s through 

(Upfm = {W'p{k)W'p{k)f'f{k) , (7.13) 



and to define 
Expanding (17.111) . 



Lf = -AUif. (7.14) 



^ j Jkfi{k)W%k,t)l^^ = {fi,e^'U'^fi) , (7.15) 

where (-, ■) denotes now the inner product in L'^(T^, dk). Using the definition of L, (17.151) 
can be written in the more symmetric form as 

(r), e^'Ulrf) = {Upv, exp[-U^'LU^h]U^fi) . (7.16) 

From the linearization one obtains 

{Lf){ki) = n dk2dk'idkA5{ki + k2 — k^ — k4)5{uJi + UJ2 — uj^ — UJ4) 

X \V{k2 - ks) + 9V{k2 - k^)\''Wl,Wl2Wl^WlS\ + h-h- U) ■ (7.17) 
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The quadratic form associated to L reads then 



TX 



4 



(/, Lf) = — dkidk2dk2,dki5{ki + k2 — k^ — k4)6{uji + iU2 — lus — ^4) 



(Td) 
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X \V{k2 - k,) + eV{k2 - k,)\^Wl,Wl2Wl,Wl,{h + f2-h- Uf ■ (7.18) 

Therefore L = L* and L > 0. Clearly LI = and Leu = 0. Thus the zero subspace is at 
least two-fold degenerate. If the H-theorem holds, see Appendix B, then the eigenvalue 
of L is exactly two-fold degenerate. In brackets, we note that momentum conservation 
is destroyed by the underlying lattice. More precisely, if one sets f{k) = k mod 1, then 
Lf 7^ 0, since the umklapp ki + k2 — k-i — ki = n,n integer vector and n 7^ 0, is permitted 
according to the momentum 5-function. 

Allowing for the interchange of limits we arrive at 

Conjecture 7.1 Under suitable conditions on uj it holds 

\imCl{\-H) = {U^fi.eM-Up^LUp^nUpfl) ■ (7.19) 

Duhamel expansion. We want to explore whether, at least in principle, the Duhamel 
expansion of Section [21 could work. The starting point is (12.281) for hq = 2. This results in 

^O'(^) = ^o'main(^) + ^oVror(i) • (7.20) 

The main term is discussed first. 

We insert the main term in the definition of Cj{t). Then the time i = input for the 
Feynman diagrams is of the generic form 

n rj n 

J2((H:,ll^ik,,a,)))p,, = ^(n«(^.'^.)W(^)L=o- (7.21) 

Therefore we have to require a slightly modified £i-clustering as 

Assumption A5. (strengthened £i-clustering). Let /3 > and /i satisfy d<5.iD . There exists 
Ao > and cq > independent ofn such that for < A < Aq and all n > 4 one has the 
bound 

E '^x,o|^(n«(^.'^.))AA(^)L=o|^^(co)"n!. (7.22) 



In addition the following limit exists uniformly in k 

. (7.23) 



l^oE |^(«(^)*«(o))/^.^(^)L=o - (u'^vYix) 



x&L 

Under Assumption A5, the input function for the Feynman diagrams has the same 
regularity as used in Section[6l We conclude that together with Assumptions Al, A3, and 
A4, the Duhamel expansion converges term to the conjectured limit (17.191 ). 
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For the error term one has symbolically 

H:i,erroAt) = f dsJ^^Jt - s) [a(s)] . (7.24) 

For the Kubo inner product we use Schwarz inequality as in (16.391 ) to conclude 

/ . {{HZ, -f^0,crror(^)))/3,A 



uiG, 



7d 



<t{J2{{HlH3,,,) / dsJ2{{n,oit-sMs)],J^lJt-s)[a{s)])),,,.a.25) 

The first factor is bounded. In the second factor we use stationary to arrive at 

/ ds y^ {{J='N,(}{s)[a\, TN,w{s)[a])) i3,x . (V.26) 

^0 „.^^, 



WdlL 



Proceeding as before, one would like to take the ^^ in the exponential to modify the state 
{■)ji,\- But this would be a high order polynomial, hence difficult to control. One may 
view (17.261) also as Feynman diagrams with singular initial conditions resulting from the 
one 5-function because of the ^^. This would require to redo the oscillatory integrals. 
At present it is not clear whether such an approach could work even in principle. 

A Appendix. Truncated correlation functions 

We follow Bratelli and Robinson IfTSl . pages 39 and 43. 

Let (■) be an even state on the CCR, resp. CAR, algebra. The moments are assumed 
to exist and odd moments vanish. We use aj as shorthand for (/j, ctj(o"j)). Then the 
moments are 

n 

(J] a,) (A.l) 

i=i 
and the fully truncated moments are denoted by 



n 



l^a.f. (A.2) 

To define them, let Z denote an index set and F a function from the non empty ordered 
subsets of 5 to the complex numbers. The truncation Ft of F is now given recursively by 

F{l) = Y,e{Vi) WFr{J), (A.3) 

where the sum is over all partitions Vi of / into ordered even subsets, Vi = { Ji, . . . , J„} 
and £(Vj) = 1 for bosons and e(Vj) = ±1 is according to whether the permutation 
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/ 1-^ (Ji, . . . , Jn) is even or odd. In our case F are the moments and Ft their full 
truncation. 

Theorem A.l (truncated correlations). Let n be even, n> A, and let n be a permutation 

of I = (1, . . . , n). Then 

n n 

{l[a,f = e{7r){l[a^^,^f, (A.4) 

i=i i=i 

where eiji) = Ifor bosons and for fermions eij:) = signvr, the sign of the permutation vr. 

Proof: Let / = {1, . . . , n} and a = {m, m + 1} C /. We set 

F{I) = (ai . . . a„) , F{I) = {a^ . . . a^+ia^ • • • a„) . (A.5) 

It holds 

F{I) = Y, <J^ I ^ J)Ft{J)F{I \ J) , (A.6) 

aeJCI 

where e{J, I \J) = 1 for bosons and £:( J, / \ J) = ±1 for fermions according to whether 
/ ^— *> {J, I --^ J) is an even or odd permutation of /. Since amam+i — ^Om+i«m = cl, 
c G C, one has 

cF{I\a) = F{I)-eF{I)= ^ e(J,/\ J)(FT(J)-0Fr(J))F(/\ J) (A.7) 

aeJCl 

= {FT{a)-eFT{a))F{I\a)+ ^ e(J, / \ J)(Fr(J) - ^Ft(J))F(/ \ J) . 

Since Fxia) — OFria) = c, we conclude 

0= J2 e{J,I\J){FT{J)-eFT{J))F{I\J). (A.8) 

For \I\ = 4, the only summand is J = / and Ft{I) = 9Ft{I) since F(0) = 1. By 
iterating (IA.8|) . one concludes the validity of (|A.4|) . □ 



B Appendix. Some properties of the spatially homoge- 
neous Boltzmann-Nordheim equation 

The purpose of this appendix is to provide a rather compressed list of the basic properties 
of the spatially homogeneous Boltzmann-Nordheim equation and to point at the relevant 
literature. It will be convenient to discuss fermions and bosons separately. 
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B.l Fermions 

The Boltzmann-Nordheim equation reads, see (15.41) . (15.51) 

d 



dt 



Wik,t)=C{W{t)){k), (B.l) 



C{W){ki) = n dk2dksdk4S{ki + k2 — k^ — k4)5{uji + uj2 — ^3 — ^4) 

\V{k2 - ks) - V{k2 - ^4)^(1^1 1^2 ^^31^4 - W1W2W3W4) , (B.2) 

W = 1 — W.lt has to be solved for initial data W such that < W^(A;) < 1 . If there would 
be a first time t such that W{ko,t) = 0, then C(W{t)){kQ) > 0, and correspondingly 
W{ki,t) = 1 implies that C{W{t)){ki) < 0. Hence the constraint < VT < 1 is 
preserved in time. Note that C(W) = in case V{k) = V{0). 
At least formally, total number and energy are conserved, 

[ dkiC{W){ki)=0, f dkiuj{ki)C{W){ki)=0. (B.3) 

According to the Fermi statistics, the entropy per volume of a quasifree state is given by 

S{W) = - f dk {W (k) logW {k)+W{k) log W{k)). (B.4) 

The entropy changes in time as 



with the entropy production 



j^S{W{t))=a{W{t)) (B.5) 



a{W) = n dkidk2dksdk45{ki + k2 — k^ — k4)6{uJi + uj2 — ujs — UJ4) 

\V{k2 - ks) - V{k2 - k4)\'^F{WiW2W3W4, W1W2W3W4) (B.6) 

and F{x, y) = {x — y) \og{x/y). Hence 

a > , (B.7) 

which is the H-theorem. a = if and only if 

W1W2 W3W4 



W1W2 W3W4 



(B.8) 



on the collision set Vc = {(^i, k2, k^, k^) E 9}'^\ki + ^2 = ^3 + ^4, coi + uj2 = UJ3 + UJ4}. 
If one introduces 

= hg{W/W) , (B.9) 
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then is a coUisional invariant in the sense that 

01 + 02 = 03 + 04 OnPe- (B.IO) 

Under some regularity on uj and for J^^ dk\(j){k) | < oo, it is proved in [l32ll that the only 
solutions to (IB.IOI ) are 

(P{k) = a + hu{k). (B.ll) 

If ly is stationary and if J^^ dk\ log (W^/(l — W))\ < oo, then W is necessarily of the 
form 

W^^^{k) = (e^(-W-^) + 1)"^ (B.12) 

for some /3, /i G M. In the limit /5 — > oo, one obtains 

iy-,^(fc) = lKfc)</i) (B.13) 

and for /? -^ — oo 

Wl^^^ik) = liujik) > fi) . (B.14) 

Although the corresponding collision invariant is not integrable, one checks directly that 
W^^^ and Wl^^^ are stationary solutions of (IB.ll ). ( IB.2I ). Presumably our list comprises 
all stationary solutions of the Boltzmann-Nordheim equation. 
Let us introduce the density, p, and energy, e, of W through 

p{W)= f dkW{k), e{W)= f dkuj{k)W{k). (B.15) 

If one picks some W, then there is a unique pair /?, /i such that p{W) = p(VF^ ), e(Vr) = 
e(M/r ), provided one also admits the values [3 = ±oo. Because of the conservation laws, 
this strongly suggests that for the solution with this W as initial datum it holds 

lim ly (t) = ly- . (B.16) 



t^oo 



On the mathematical side, besides the study of collision invariants, the main focus 
so far is the existence and uniqueness of solutions. We refer to the review [|33l . These 
authors consider the momentum space M'^ and dispersion uj(k) = k"^, see Remark 5.3. To 
be concise we quote one result from Dobeault [ j34l . 

Theorem B.l Let d > 2 and uj{k) = k"^. Let V be radial with J^^ dk\k\\V {k)\'^ < oo. 
For the initial datum, W, we assume W G L°°{W^) and < py < 1. Then the integrated 
version of the Boltzmann-Nordheim equation, 

W{t) = W+ f dsC{W{s)) , (B.17) 

Jo 

has a unique solution in C([0, oo), L°°{W^)), which satisfies < W{t) < Ifor all t > 0. 
IfpiW) < oo, e{W) < oo, and S{W) < oo, then it holds, for all t >0, 



dkW{k,t) = p{W), / dkkW{k,t)= dkkW{k), 
f dkk'^W{k,t) = e{W), S{W{t)) = S{W) + f dsa{W{s)) . (B.18) 

JRd Jo 
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Apparently, the case of interest in our context has never been investigated. One diffi- 
culty results from the constraint due to the energy (5-function. On the other hand, since T*^ 
is compact, finite number, energy, and entropy holds automatically. Under the Assump- 
tion A2 we expect Theorem B. 1 still to be valid. 

B.2 Bosons 

The Boltzmann-Nordheim equation reads, see (|5.4I) . (|5.5I) . 

yV{k,t)=C{W{t)){k), (B.19) 



CiyV){k) = TT dk2dk^dk4^5{ki + k2 — k^ — k4)6{iUi + u!2 — ujs — UJ4) 

\V{k2 - k3) + V{k2 - ki)\^{WiW2W3Wi - W1W2W3W4) , (B.20) 

W = 1 + W. It has to be solved with initial data W such that W > 0, f dkW{k) < 00. 
Positivity is preserved in the course of time and number and energy conservation holds as 
in (131). 

According to the Bose statistics, the entropy per volume of a quasifree state is given 
by 

S{W) = - I dk{W{k)\ogW{k)-W{k)\ogW{k)). (B.21) 

Then ^3 to (IBUI) still hold, provided W = l + W widm ^M \V{k2 - k^) - V{k2 - 
ki)\'^ is substituted by \V{k2 - k^) + V{k2 - ^4)^. 

If W is stationary and j^^ dk\ log W^/(l + W)\ < 00, then W is necessarily of the 
form 

Wl^{k) = (e^(-W-^) - 1)"\ (B.22) 

where, in order to have Wt„ > 0, one is restricted to 

/3 > , CJmin > yU and P <0, fl > CUmax ■ (B.23) 

The density and energy of a state W is still given by (IB. 151 ). Clearly, the range of this 
map is D = {{p, e)|0 < p < cxo, pcUmin < e < pcjmax, }• However for rf > 3 and under 
our assumptions on u, the map (/5, p) -^ (p(W^/3^^); ^(PV^^)) has a range Dnov which is 
a proper subset of D. Therefore, for an initial W with {p{W), e{W)) E D \ Dnor the 
solution of the Boltzmann-Nordheim equation has no limiting stationary state. Physically 
the excess mass condenses to a superfluid component, which is reflected by a (5-peak 
at /c = in the Wigner function. Based on self-similar solutions it is argued that the 
condensate is nucleated at some finite time. We refer to the review ll35l . see also |[36l[37l . 

On the mathematical side, the cubic nonlinearity of the Boltzmann-Nordheim equa- 
tion poses severe difficulties. Only the case of momentum space M'^, dispersion u{k) = 
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k"^ /2 has been studied in detail and even then only for isotropic solutions which makes 
W{k^ t) to depend on |A;| only. We use e as energy variable, e = k'^/2, and set W{k, t) = 
f{e,t). In case y (A;) = V{0) = l/(2(27r)^) the Boltzmann-Nordheim equation for the 
energy distribution reads 

^/(£i,t)=Cr(/(t))(5i) (B.24) 

with the collision operator 

Cr(/)(£i) = / de3dei^= min { v^, v^, y/e^, V^ }(/i/2/3/4 - /1/2/3/4) , 

(B.25) 
see [|36l . Here 62 = 6^ + e^ — ei, fj = f{ej), j = 1, . . . , 4, / = 1 + /. Lu [|38l considers 
Equation (IB .241) in the weak form 

-J ^{e)f{e,t)^ede = j ^{e)Cr{f{t)){e)^ede (B.26) 

with test functions ip E Cl(K+) and proves that (IB.26I) remains meaningful even if 
f{e,t)y/ede is a positive measure on M+. We denote such a measure by f{de,t). The 
possible roughness of /{de, t) is balanced by rewriting the right hand side of (IB.26I) in 
such a way that the coUisional difference ^{ei) + ^{£2) — ^(^3) — ^i^i) appears. For 
t = we assume finite mass and energy, i.e., J^ f{de) < 00 and J^ ef{de) < 00. 
Then (IB. 261) has a solution which conserves mass and energy. The stationary measures of 
(IB. 261) are necessarily of the form 

fp,^{de) = (e^(^-^) - ly'V^de + n,,J{e)de (B.27) 

with either /i < and condensate density ricon = or yU = and ricon > 0. The critical 
line dividing the normal fluid from the condensate is given by e(p) = (p/ po)~^/^eo with 
po = r deV^ie^^ - l)-\ eo = J,°° de^ee{eP^ - l)-\ 

For the particular case (IB. 241) one has some rigorous information on the concentration. 
Let the initial data be given by a density, i.e., f{de) = fac{e)^de, such that 



P = pU) 



/•OO /"CO 

/ f.,c{e)\fede < cx) , e = e(/) = / ef^{e)y^de < cx) . (B.28) 
Jo Jo 



In general, the solution to (|B.26I) is then a measure, f{de,t), which can be uniquely 
decomposed into an absolutely continuous and a singular part, 

fide,t) = U,{e,t)^ede + f,{de,t) . (B.29) 

Now let fi3^^ be the equilibrium distribution corresponding to p(/), e(/), where /i = in 
case neon > 0. Lu 1,39,1 proves that, for p < Pc, 

/•OO /"OO 

lim / /3(ci5,t) = 0, lim / \U,{e,t)-fp^^{e)\^ede = Q. (B.30) 
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On the other hand for p > Pc one has 

POO 

Urn / \f^{e,t)-fp,oie)\eVEde = (B.31) 



and 



lim / f{de,t)=p-p,, (B.32) 



t^oo l{e\0<ei<r{t)} 



where r{t) is the integral in (|B.31I) . Thus for large times the solution develops a 5-peak at 
e = with weight ricon- 

Whether such a concentration appears after some finite time or only asymptotically, 
as t ^ oo, is left unanswered by the theorem. 
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